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ABSTRACT 

We present unambiguous evidence for a parsec scale wind in the Broad- 
Line Radio Galaxy (BLRG) 3C 382, the first radio-loud AGN, with = 
l°Sio(/5GHz//44oo) > 1, whereby an outflow has been measured with X-ray grat- 
ing spectroscopy. A 118ks Chandra grating (HETG) observation of 3C382 has 
revealed the presence of several high ionization absorption lines in the soft X-ray 
band, from Fe, Ne, Mg and Si. The absorption lines are blue-shifted with re- 
spect to the systemic velocity of 3C382 by —840 ±60kms _1 and are resolved by 
Chandra with a velocity width of a = 340 ± TOkms^ 1 . The outflow appears to 
originate from a single zone of gas of column density Ah = 1.3 x 10 21 cm -2 and 
ionization parameter log(£/ergcms -1 ) = 2.45. From the above measurements 
we calculate that the outflow is observed on parsec scales, within the likely range 
from 10 — 1000 pc, i.e., consistent with an origin in the Narrow Line Region. 



Subject headings: Galaxies: active — galaxies: individual (3C382) — X-rays: 
galaxies 
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1. Introduction 
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tion and emission features at soft (0.4-2 keV) and hard (6-8 keV) X-ray energies, coincident 
with ionized O, N, Ne, Mg, Si and Fe lines blueshifted in the observer's rest-frame. The in- 
ferred velocities of the win ds are typically in the range 100 — 1000 km s~\ but can be as high 



as ~ 0.1c in some sources ( Chartas et al.ll2002l . 12003 



Pounds et al. 2003; Reeves et al. 2003 



20081 ). It is also possible that the energy budget of the outflows of so me AGN can approac h 



a significant fraction of the bolometric or even Eddington luminosity (IKing fe Pounds 



2003|). 



In stark contrast, the X-ray evidence for nuclear outflows is very scarce in Broad-Line 
Radio Galaxies (BLRGs) and in radio-loud AGN generally. Previously, th e radio-loud quasa r 
4C +74.26 showed weak absorption features at ~ 1 keV with ASCA ( IBallantynd 120051 ). 
while the BLRG Arp 102B showed neutral X-ray absorption with ASCA and a UV outflow 



of a few hu ndred kms 1 (lEracleous et al. 
and 3C 33 ( Sambruna et al. 20071 : lEvans et ah 



2003). Furthermore two radio galaxies, 3C 445 



20061 ) also exhibit soft X-ray emission lines 



below 2 keV, which could originate from spatially extended material (in 3C 33; iTorresi et al. 



(2009a 


_ ~ — , , .. „„„ _ 
)). Disk wind 


3 are 


(Blandford & Pavne 


1982) 



To differentiate between radio-loud and radio-quiet AGN, here we adopt the radio- 
loudness parameter R L = log 10 (/ 5 GHz//44oo), where f^aw* is the cor e 5 GHz radio flux and 
/4400 is the flux at 4400A, both in units of mJy ( Kellerman et al.lll989). Generally, Ri, > 1 fo r 
radio-loud AGN Jwilkes fe Elvidll987h . while for 3C382, R L = 1.9 Lawson fe Turner!E997l ). 
If the extended radio emission from 3C382 is also included, then may be considerably 
higher. 

In this Letter we present direct evidence for outflowing gas from the nucleus of the 
nearby {z = 0.05787), bright BLRG 3C 382. A r e-analysis of our 1 18 ks Chandra HETG 
(High Energy Transmission Grating) observations (IGliozzi et al.ll2007l ) revealed several blue- 
shifted absorption lines between 0.7 — 2.0 keV which suggest the presence of a large-scale 
(10 — 1000 pc) outflow in this source with a velocity of 800 kms -1 . The organization of this 
Letter is as follows. In § 2 we describe the Chandra data reduction and analysis; in § 3 
the results of the spectral analysis; Discussion and Conclusions follow in § 4. Throughout 
this paper, a concor dance cosmology with H = 71 km s -1 Mpc -1 , ^a=0.73, and f2 m =0.27 
( jSpergel et al.l 120031 ) is adopted. Errors are quoted to 90% confidence for 1 parameter of 
interest (i.e. Ax 2 or AC = 2.71). 
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2. The Chandra HETG Data 

Chandra observed 3C382 with the HETG for a net exposure of 118 ks between 27-30 
November 2005. The ±1 order spectra were summed for the MEG (Medium Energy Grating) 
and HEG (High Energy Grating) respectively, along with their response files. The summed 
first order count rates for the MEG and HEG are 0.867 counts s _1 and 0.379 counts s" 1 respec- 
tively, while the MEG data were fitted between 0.5-7.0 keV and the HEG from 1.0-9.0 keV. 



3. The Warm Absorber in 3C 382 



The Chandra HETG data were first analyzed by iGliozzi et al.l (120071 ). who focused on 
the continuum and its variability. Our results are in agreement with theirs. Specifically, 
the MEG and HEG data were fitted by an absorbed power law with photon index T = 
■66±0.01plusab 



ackbody with kT = 92 ±6 eV to parameterize the soft excess below 1 keV 
Gli ozzi et al.ll2007t). absor bed by a Galactic line of sight column of A^H,Gai = 7.0 x 10 20 cm~ 2 
( iDickey fc Lockmanl Il990l ). Figure 1 shows the broad-band HETG spectrum fitted with 



an absorbed power-law only, to illustrate that a soft excess is clearly present below 1 keV. 
The 0.5-9 keV band flux is 6.4 x 10 -11 erg cm -2 s -1 . Even upon adding a blackbody to 
parameterize the soft excess, the fit is still formally unacceptable (x 2 /dof = 632/477, null 
probability 2.4 x 10~ 6 ) as there are clear residuals around 1 keV that indicate the presence 
of a warm absorber. 

To analyse the warm absorber in detail, the HEG and MEG spectra were binned more 
finely to sample the resolution of the detector, at approximately HWHM the sp ectral resolu - 
tion (e.g. AA = 0.01 A bins for the MEG). For the fits, the C-statistic was used dCashlll979h . 
as there are fewer than 20 counts per resolution bin. The absorption lines were modelled 
with Gaussian profiles and the continuum model was adopted from above. Table 1 lists the 
detected lines with their observed and inferred properties, and their significance as per the 
C-statistic. Figure 2 shows the portions of the HETG spectrum containing the strongest 
lines, with the model overlaid. 

The seven absorption lines in Table 1 and Figure 2 are all detected at high confidence 
(corresponding to AC > 18, or > 99.9% confidence for 2 parameters of interest). The lines 
likely arise from the Is — 2p transitions of Neix, Nex, Sixm, and Sixiv and the 2p — 3d 
lines of Fexix-xxi. The two statistically weaker Is — 2p lines of Mgxi and Mgxn may also 
be present, which have outflow velocities consistent with the other lines. 

Initially we assume that the lines have the same velocity width within the errors. The 
velocity width of the absorption lines is then o = 340±70kms~ 1 (or 780±160kms~ 1 FWHM) 
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and the lines are clearly resolved. Even at 99% confidence (AC = 9.2 for 2 parameters), the 
velocity width is constrained to a = 340±140 km s _1 . Upon allowing the velocity width of the 
individual lines to vary, then they are constrained to lie within the range a = 250— 500 km s -1 
as shown in Table 1. The mean outflow velocity is —810 km s -1 . The overall fit statistic is 
C = 3804 for 3811 bins. 



We used the photoionization code xstar ( jKallman et al.ll2004l ) to derive the parameters 



of the absorber, assuming the baseline conti nuum described above, in cluding the soft excess. 
Solar abundances are assumed throughout ( Grevesse fc Sauval 1998 ). An important input 



parameter is the turbulent velocity, which can effect the absorption line equivalent widths and 
hence the derived column density. We experimented with two different values of the turbulent 
velocity chosen to represent two likely extremes: (i) a lowest value of t> tU rb = lOOkms -1 
and (ii) w tU rb = 300 km s^ 1 , the latter being consistent with the measured width of the 
absorption lines. The fitted continuum parameters are T = 1.68±0.02 and for the blackbody, 
kT = 110 ± 8eV. For case (i), then iVn = (3.2 ± 0.6) x 10 21 cm -2 , the ionization parameter 
is log^0= 2.45t[J;o8 an d the outflow velocity is v out = — 810^° kms -1 . The fit statistic is 
C/bins = 3795/3811. For case (ii), then N H = (1.30 ± 0.25) x 10 21 cm" 2 , log£ = 2.45lg;8? 
and the outflow velocity f out = — 840lgQ kms -1 . The fit statistic is C/bins = 3783/3811. If 
the warm absorber is not included in the model, then the fit statistic is substantially worse 
by AC = 220 (compared to model (ii)). Only a single outflowing layer of gas is required to 
model the warm absorber. 

The higher turbulence velocity model is statistically preferred and is consistent with 
the measured 340kms _1 widths of the lines. Either model yields an outflow velocity of 
— 800 km s -1 within a statistical error of < 10%. Hereafter we adopt the parameters from 
model (ii), as the turbulent velocity is consistent with widths of the individual absorption 
lines. However in neither model is the fit ted column density of the absorber as high as the 



value of A^h ~ 3 x 10 22 cm 2 reported by iTorresi et al.l (j2009bl ) from an analysis of a short 



34.5 ks XMM-Newton/RGS observation of 3C382 on April 28, 2008. If the column density 
of the warm absorber is fixed to the value of A^h = 3 x 10 22 cm~ 2 in the HETG spectrum, 
then the fit statistic is considerably worse (C/bins = 8390/3811). 

As a consistency check, we analyzed the archival RGS data of 3C382 with model (ii) 
and the same continuum form as above. We found that the column density and ionization 
parameter were degenerate with each other, given the short exposure of the RGS spectrum. 
Thus for the RGS, N u = 1A+H x 10 22 cm" 2 , log£ = 3.4 ± 1.0 and the outflow velocity 
v out = -1 200^00 km s" 1 . The fit statistic improves only by AC = 25 to C/bins = 826/820 



The units of £ are erg cm s 
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upon adding the absorber. Thus within the larger errors, the parameters are consistent 
with those obtained from the HETG. The individual lines detecte d in the Chandra HE TG 
observation were also compared to the absorption lines claimed by iTorresi et al.l (l2009bl ) on 
the basis of the RGS data. With the higher signal to noise ratio of the HETG spectrum 
compared to the RG S spectrum, we only confirm the detection of one line reported by 
Torresi et al.l ( 2009bl ). the Fexx line at a rest-frame energy of 1025 eV@ 



4. Discussion and Conclusions 

The Chandra HETG spectra have revealed an ionized outflow in the BLRG 3C382. 
The outflow parameters are well determined, with Ah = (1.30 ± 0.25) x 10 21 cm -2 , log£ = 
2.45lo'gy an d v out = —840^5Q km s" 1 , while the absorption line widths are resolved with 
a = 340 ±70 km s" 1 . 

To characterize the outflow, we define the (unabsorbed) ionizing luminosity L ion , which 
in xstar is defined from 1 — lOOORydberg. This depends on the continuum model fitted 
to the data and it is important to take into account the soft excess. Using the best fit 
powerlaw plus blackbody continuum, then L ion = 1.2 x 10 45 ergs _1 . We note that if a 
different continuum form is used to parameterize the soft excess, e.g. a broken power-law, 
then the ionizing luminosity can be a factor of ~ 2 higher; keeping this caveat in mind we 
adopt the more conservative lower luminosity value of 1.2 x 10 45 ergs _1 



4.1. The Location of the Absorber 

The upper-bound to the wind radius (R on t) is determined by geometrical constraints, 
i.e. if the thickness of the absorber AR/R << 1 (valid for a thin shell) and as Ah = nAR, 
where n is the electron number density. As the ionization parameter of the absorber is defined 
as £ = L ion /nR 2 , then by substitution R out « L ion /Nii^ = 3.3 x 10 21 cm (or << lkpcjE 
The lower bound is formed by the escape velocity, i.e. for the gas to escape the system as an 
outflow then R esc > c 2 /v 2 R s , where R s = 2GM/c 2 is the black hole Schwar zschild radius and 
v = 800 km s" 1 . If for 3C382, M = 1 x 10 9 M (with a 40% uncertainty, see iMarchesini et al. 



(120041 ) ). then R esc > 1.4 x 10 5 R S > 4.2 x 10 19 cm > 13pc. In other words the location of the 



2 The line reported at 1.356 keV by lTorresi et al.1 (|2009bl ) may also be associated with the Is — 2p line of 
Mgxi, as noted in Table 1. 

3 Note the same radius is obtained by integrating down the line of a sight of a homogeneous radial outflow. 



- 6- 



soft X-ray outflow is likely bounded between approximately 10 pc and 1 kpc. 



Furthermore the outflow velocity and absorption line FWHM of ~ 800 km s -1 is similar 
to the width of the narrow optical forbidden lines of [O ill] , [O i] and [ Si ill , which for 3C 382 
l ie in the range 400 — 600 km s -1 , as measured fro m the spectra of lEracleous fc Halpern 
( 11994J ). In particular the [O ill] emission line from the lEracleous fc Halpern! (119941 ) spectrum 
appears to have an asymmetric profile, with the blue-wing extending to ~ — 1070kms _1 
from the line centroid, which appears ~ 370 km s -1 broader than the red-wing. 

Although this is not evidence for observing outflowing gas through the direct line of 
sight, the [Olll] emission may suggest we are viewing outflowi ng gas along a different sight - 
line. Indeed this effect is also seen in other radio galaxies (IGelderman fc Whittle! Il994l ). 
The coincidence between the soft X-ray absorbing ga s and extended [O ill] emission has also 
been noted in the radio-quiet quasar, MR 2251-178 (IKaspi et al.l 120041 ) . Thus the origin of 
the X-ray absorption appears consistent with the optical Narrow Line Region (NLR). The 
association between the soft X-ray absorption in 3C 382 and any rest-frame UV absorption 
could be tested with future simultaneous Chandra and HST observations. 



The density of the outflow is then n = L\ on /^R 2 . For the parameters above, then 
0.4 — 2400 cm -3 . While loosely c onstrained, this is consistent with typical expected 



n = U.4 — z4UUcm . While loosely c onstrained, this is consistent with typical expected 
NLR densities of ~ 10 3 cm -3 in AGN fcoskilll978h . The mass outflow rate for a uniform 



spherical flow is M out = 4irnR 2 m p v ont , where nR 2 = L- lon /^ and m p is the proton mass. 



Hence for the measured warm absorber parameters for 3C382, M out 
M out = 100 Moyr" 1 . 



7.2 x 10 27 gs 



or 



However the low rate of appearance of such absorbers amongst BLRG's or radio-loud 
AGN generally might suggest that the solid angle subtended by individual absorbing fila- 
ments is much smaller than 47r steradians. Therefore, the mass outflow rate could be con- 
siderably smaller than the above estimate. Subsequently the kinetic power of the soft X-ray 
outflow (for f out = — 800 km s" 1 ) is then E = l/2M OVLt v 2 ut < 2 x 10 43 ergs _1 , which ener- 
getically is a fairly insignificant 2% of the ioniz ing luminosity and is unlikely to contribute 
significantly towards AGN feedback (lKing!l2003l ). 
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Fig. 1— Chandra HETG MEG (black) and HEG (red) spectra of 3C382, binned coarsely 
at four times the resolution of the gratings, in order to show the broad-band continuum 
spectrum. The data is plotted as a ratio against an absorbed power-law, of photon index 
T = 1.66. A clear excess of counts is present in the soft X-ray band below 1 keV, while 
significant residuals below unity are also present between 0.7-1.0 keV, indicating that a 
warm absorber is present. 
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Fig. 2. — Fluxed HETG spectra showing the comparison between the data binned at HWHM 
resolution and the best-fit absorption line model (solid line) described in the text. Several 
absorption lines are present in the HETG spectrum, as labelled in the above figure and listed 
in Table 1. 
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Table 1. Summary of HETG absorption line parameters. 



E(obs) a E(rest) 6 ID C EW rf v e out a? AC<> 



868.5 918.8±J;^ Fexix 2p - 3d (917.0) -1.4 ±0.5 590±3gg 509^™ 24 
874.1 924.7±g;? Neix Is - 2p (922.0) -1.9±^;| 875±||^ 493j^2 35 
915.8 968.8±g;3 Fe xx 2p - 3d (967.3) -1.6 ±0.5 460±^ 516j^g 21 

956.6 1012.01^ Fe xxi 2p - 3d (1009.0) — 1.4±°;| 890 ± 210 296+^° 19 
969.0 1025. l±oi Nex Is - 2p (1021.5) — 2.1±g;| 1050 ± 120 289±^| 48 
1767.7 1870.0l^j Si XIII Is - 2p (1865.0) — 1.6±g;| 800±^ 255+^ 26 
1901.1 2011.ll?j Si xiv Is - 2p (2004.4) -2.2 ±0.7 1000±fg 284+^ 32 
1281.4 1355.6±i;g Mgxi Is - 2p (1352.2) -1.0 ±0.4 750 ± 220 340 h 12 
1395.4 1476.21^ Mgxn Is - 2p (1472.2) -0.8 ± 0.4 810 ± 203 340 h 9 



a Observed energy of absorption line in eV. 

b Energy of absorption line in rest-frame of 3C 382, in units eV. 

c ine identification and lab frame energy in eV in parenthesis. Atomic data are from 



Equivalent width, units eV. 
e Outflow velocity of absorption line, in units kms -1 . 
f lcr velocity width of absorption line, in units kms -1 . 
g Improvement in C-statistic, upon adding line to model. 
h Parameter is fixed in the model. 



